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Abstract. The processes of radiation damping and decoherence in Quantum Elec- 
trodynamics are studied from an open system's point of view. Employing functional 
techniques of field theory, the degrees of freedom of the radiation field are eliminated 
to obtain the influence phase functional which describes the reduced dynamics of 
the matter variables. The general theory is applied to the dynamics of a single 
electron in the radiation field. From a study of the wave packet dynamics a quanti- 
tative measure for the degree of decoherence, the decoherence function, is deduced. 
The latter is shown to describe the emergence of decoherence through the emission 
of bremsstrahlung caused by the relative motion of interfering wave packets. It is 
argued that this mechanism is the most fundamental process in Quantum Elec- 
trodynamics leading to the destruction of coherence, since it dominates for short 
times and because it is at work even in the electromagnetic field vacuum at zero 
temperature. It turns out that decoherence trough bremsstrahlung is very small for 
single electrons but extremely large for superpositions of many-particle states. 



1 Introduction 

Decoherence may be defined as the (partial) destruction of quantum coher- 
ence through the interaction of a quantum mechanical system with its sur- 
roundings. In the theoretical analysis decoherence can be studied with the 
help of simple microscopic models which describe, for example, the interaction 
of a quantum mechanical system with a collection of an infinite number of 
harmonic oscillators, representing the environmental degrees of freedom 0,0. 
In an open system's approach to decoherence one derives dynamic equations 
for the reduced density matrix || which yields the state of the system of 
interest as it is obtained from an average over the degrees of freedom of the 
environment and the resulting loss of information on the entangled state of 
the combined total system. The strong suppression of coherence can then be 
explained by showing that the reduced density matrix equation leads to an 
extremely rapid transitions of a coherent superposition to an incoherent sta- 
tistical mixture [Ep|. For certain superpositions the associated decoherence 
time scale is often found to be smaller than the corresponding relaxation or 
damping time by many orders of magnitude. This is a signature for the fun- 
damental distinction between the notions of decoherence and of dissipation. 
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A series of interesting experimental investigations of decoherence have been 
performed as, for example, experiments on Schrodinger cat states of a cavity 
field mode [|| and on single trapped ions in a controllable environment . 

If one considers the coherence of charged matter, it is the electromagnetic 
field which plays the role of the environment. It is the purpose of this paper to 
study the emergence of decoherence processes in Quantum Electrodynamics 
(QED) from an open system's point of view, that is by an elimination of 
the degrees of freedom of the radiation field. An appropriate technique to 
achieve this goal is the use of functional methods from field theory. In section 
2 we combine these methods with a super-operator approach to derive an 
exact, relativistic representation for the reduced density matrix of the matter 
degrees of freedom. This representation involves an influence phase functional 
that completely describes the influence of the electromagnetic radiation field 
on the matter dynamics. The influence phase functional may be viewed as 
a super-operator representation of the Feynman- Vernon influence phase [jjj 
which is usually obtained with the help of path integral techniques. 

In section 3 we treat the problem of a single electron in the radiation 
field within the non-relativistic approximation. Starting from the influence 
phase functional, we formulate the reduced electron motion in terms of a 
path integral which involves an effective action functional. The corresponding 
classical equations of motion are demonstrated to yield the Abraham-Lorentz 
equation describing the radiation damping of the electron motion. In addition, 
the influence phase is shown to lead to a decoherence function which provides 
a measure for the degree of decoherence. 

The general theory will be illustrated with the help of two examples, 
namely a free electron (section 4) and an electron moving in a harmonic po- 
tential (section 5) . For both cases an analytical expression for the decoherence 
function is found, which describes how the radiation field affects the electron 
coherence. 

We shall use the obtained expressions to investigate in detail the time- 
evolution of Gaussian wave packets. We study the influence of the radia- 
tion field on the interference pattern which results from the collision of two 
moving wave packets of a coherent superposition. It turns out that the ba- 
sic mechanism leading to the decoherence of matter waves is the emission 
of bremsstrahlung through the moving wave packets. The resultant picture 
of decoherence is shown to yield expressions for the decoherence time and 
length scales which differ substantially from the conventional estimates de- 
rived from the prominent Caldeira-Leggett master equation. In particular, it 
will be shown that a superposition of two wave packets with zero velocity 
does not decohere and, thus, the usual picture of decoherence as a decay of 
the off-diagonal peaks in the corresponding density matrix does not apply to 
decoherence through bremsstrahlung. 

We investigate in section 6 the possibility of the destruction of coherence 
of the superposition of many-particle states. It will be argued that, while the 
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decoherence effect is small for single electrons at non-relativistic speed, it is 
drastically amplified for certain superpositions of many-particle states. 
Finally, we draw our conclusions in section 7. 

2 Reduced Density Matrix of the Matter Degrees of 
Freedom 

Our aim is to eliminate the variables of the electromagnetic radiation field 
to obtain an exact representation for the reduced density matrix p m of the 
matter degrees of freedom. The starting point will be the following formal 
equation which relates the density matrix p m (tf) of the matter at some final 
time tf to the density matrix p(U) of the combined matter-field system at 
some initial time U, 

p m (t / )=tr / |T^exp J ' tfxC{x) pfe) j . (1) 

The Liouville super-operator C(x) is defined as 

C(x)p = -i[H(x),p], (2) 

where TL{x) denotes the Hamiltonian density. Space-time coordinates are 
written as x^ = (a; , a;) = (t, x), where the speed of light c is set equal 
to 1. All fields are taken to be in the interaction picture and TV- indicates 
the chronological time-ordering of the interaction picture fields, while try de- 
notes the trace over the variables of the radiation field. Setting % = c = 1 we 
shall use here Heaviside-Lorentz units such that the fine structure constant 
is given by 

To be specific we choose the Coulomb gauge in the following which means 
that the Hamiltonian density takes the form ]§lfl|lO| 

H(x)=H c (x)+H tT (x). (4) 

Here, 

H tT {x) = f{x)A^x) (5) 

represents the density of the interaction of the matter current density j^{x) 
with the transversal radiation field, 

A"(x) = (0, A(x)), V • A{x) = 0, (6) 



and 



(7) 



4 



Heinz-Peter Breuer and Francesco Petruccione 



is the Coulomb energy density such that 
Hc(x°) = \(d?xjd 



3 j°(x ,x)j (x°,y) 

"i 

Att\x — y\ 



y- 



(8) 



is the instantaneous Coulomb energy. Note that we use here the convention 
that the electron charge e is included in the current density j^(x) of the 
matter. 

Our first step is a decomposition of chronological time-ordering operator 
T<_ into a time-ordering operator T-L for the matter current and a time- 
ordering operator for the electromagnetic field, 



T = T A 
This enables one to write Eq. (Q) as 



Pm(t f )=Tl try T^exp 



d 4 x (C c (x) + C tI (x)) 



p(u) 



(9) 



(10) 



where we have introduced the Liouville super-operators for the densities of 
the Coulomb field and of the transversal field, 



£ c {x)p = -i[Hc(x),p], Ctr(x)p = -\[f{x)A tl {x),p}. 



(11) 



The currents j 71 commute under the time-ordering T£_. We may therefore 
treat them formally as commuting e-number fields under the time-ordering 
symbol. Since the super-operator Cc{x) only contains matter variables, the 
corresponding contribution can be pulled out of the trace. Hence, we have 



Pm(tf) = T'i_ exp 



d A xCc(x) 



tr/ <^ T^exp 



d A xCtx{x) 



p(ti 



(12) 

We now proceed by eliminating the time-ordering of the A-fields. With 
the help of the Wick-theorem we get 



T^exp 



exp 



d 4 x£ tr (x) 



d 4 X dV[Ar(x),£tr(x')Mi™0 



(13) 







exp 


J d 4 x£ tr (x) 



In order to determine the commutator of the Liouville super-operators we 
invoke the Jacobi identity which yields for an arbitrary test density p, 

[£tr{x),£ tI (x')]p = £ tr (x)£tr{x')p- £ tr (x')£ tl (x)p 

= -[H t r(x), [H tI (x'),p}] + [Htr(x'), [H t r(x),p]} 

= -[[Htr(x),H tI (x%p]. (14) 
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The commutator of the transversal energy densities may be simplified to read 

[Ht T (x),n tI (x')}=j^x)r(x')[A^x),A v (x% (15) 

since the contribution involving the commutator of the currents vanishes by 
virtue of the time-ordering operator T£_. Thus, it follows from Eqs. Jl4| ) and 
T5|) that the commutator of the Liouville super-operators may be written as 



[C tI (x), Ar (x')]p = - [A,, (x) , A v (x')} [f(x)f (x'), p] . 



(16) 



It is useful to introduce current super-operators J+{x) and J_(x) by means 
of 

Jl{x)p = f{x) P , Jt{x)p = pf{x). (17) 

Thus, J+(x) is defined to be the current density acting from the left, while 
J-{x) acts from the right on an arbitrary density. With the help of these 
definitions we may write the commutator of the Liouville super-operators as 

[C tI (x) 7 Ctr(x')} = - [A^x),A v {x')]J»{x)r + {x') 
+ [A^x),A v {x')]Jt{x)r{x'). 



Inserting this result into Eq. (|13|), we can write Eq. ( |12| ) as 



Pm{tf) = Tl_ exp 



d 4 xCc(x) 



d 4 x 



d 4 x'9(t - t')[A^(x), A„(x')}J>t(x)Jl(x') 



d 4 x / d 4 x'9{t - t'^A^x), A,{x / )}J f l(x)r(x') 



•try < exp 



d 4 a:£ tr (2;) 



(18) 



This is an exact formal representation for the reduced density matrix of the 
matter variables. Note that the time-ordering of the radiation degrees of 
freedom has been removed and that they enter Eq. ( |l8| ) only through the 
functional 



W[J +1 J_] = tr/ < exp 



d 4 xC t r(x)p(ti) 



ti 



(19) 



since the commutator of the A-fields is a c-number function. 



3 The Influence Phase Functional of QED 



The functional (JTy) involves an average over the field variables with respect 
to the initial state p{U) of the combined matter-field system. It therefore 
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contains all correlations in the initial state of the total system. Here, we 
are interested in the destruction of coherence. Our central goal is thus to 
investigate how correlations are built up through the interaction between 
matter and radiation field. We therefore consider now an initial state of low 
entropy which is given by a product state of the form 

p(U) = Pm{U)® Pf, (20) 

where p m {ti) is the density matrix of the matter at the initial time and 
the density matrix of the radiation field describes an equilibrium state at 
temperature T, 

p f = ±exp(-PH f ). (21) 
z f 

Here, Hf denotes the Hamiltonian of the free radiation field and the quantity 
Zf = tr/[exp(— pHf)\ is the partition function with (3 = 1/ksT. In the 
following we shall denote by 

{0)f = tx f {Opf} (22) 
the average of some quantity O with respect to the thermal equilibrium state 

©' ^ 

The influence of the special choice (|20|) for the initial condition can be 

eliminated by pushing t, — > — oo and by switching on the interaction adiabat- 
ically. This is the usual procedure used in Quantum Field Theory in order to 
define asymptotic states and the S'-matrix. The matter and the field variables 
are then described as in-fields, obeying free field equations with renormalized 
mass. These fields generate physical one-particle states from the interacting 
ground state. 

For an arbitrary initial condition p(ti) the functional PF[J+, J_] can be 
determined, for example, by means of a cumulant expansion. Since the initial 
state ( p0[ ) is Gaussian with respect to the field variables and since the Liouvillc 
super-operator Ctr(%) is linear in the radiation field, the cumulant expansion 
terminates after the second order term. In addition, a linear term does not 
appear in the expansion because of (A^x)) f — 0. Thus we immediately 
obtain 



W[J+, J_] = exp 



d 4 x / dV(£ tr (z)£ tr 0E')>/ 



Pm(ti). (23) 



Inserting the definition for the Liouville super-operator Ct r (x) into the expo- 
nent of this expression one finds after some algebra, 

1 f l f ftf 

d x d x'(C tT (x)C t y;(x')}fp m 



2Jt 



' 1 ' d 4 X I f d A x'trf{[Htr(x),[H tI (x'),p m ®Pf]]} 



2 



ti Jti 
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if 



d 4 x I dV[(A„(^V.W>/^W^(^) 

+ {A^x)A v {x')) f r{x).r{x') 
-(A„(x')A^x)) f J>t(x)r(x') 

-(A^x)A„(x')) f r(x)r(x')] Pn 



On using this result together with Eq. (P3|), Eq. ( |18| ) can be cast into the 
form, 



,(*/) = Tl( 



exp 



d 4 xCc(x) 



(24) 



d 4 x / d 4 x'{-{e{t-t')[A^x),A v {x')\ 

+(A v {x')A,{x)) f )r + {x)r + {x') 

+ (6(t-t')[A^x),A l/ (x')} 

-(A,(x)A„(x<)) f )r(x)r(x') 
+ {A y {x')A,{x)) f Jl{x).r{x') 

+(A^x)A v (x')) f J^(x)Jl(x')}]) Pm (U ) . 



At this stage it is useful to introduce a new notation for the correlation 
functions of the electromagnetic field, namely the Feynman propagator and 
its complex conjugated (T_^ denotes the anti-chronological time-ordering), 

W F {x-x')^ = {T^{A„{x)A u {x'))) f 

= 0{t - t')[A^{x),A v {x')] + (A v (x')A^x)) f7 
iD* F (x - x 1 )^ = -{T^{A^x)A u {x'))) f 

= 9(t - t')[A^x), A v (x')} - (A„(x)A v (x')) f , (25) 

as well as the two-point correlation functions 



D + {x-x')^ v = {A^{x)A v (x')) f , 
D_(x - ee (A u {x')A^{x)) f . 



(26) 



As is easily verified these functions are related through 

- \D F {x - x')^ + iD* F (x - £%„ + D+(x - x')^ + D-{x - x')^ = 0. (27) 

With the help of this notation the density matrix of the matter can now be 
written as follows, 



Pm(tf) = Tl_ exp 



ftf 

/ d 4 xC c (x) 



(28) 
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+\ |''d 4 i J*' d 4 x' {~iD F (x-x')^J»(x)Jl(x') 

+iD* F {x-x')^J»(x)r{x') 

+D_(x-x') tlv J»(x)J v _{x') 

+D + (x - x')^Jt(x)JHx')}]) Pm (ti). 

This equation provides an exact representation for the matter density matrix 
which takes on the desired form: It involves the electromagnetic field variables 
only through the various two-point correlation functions introduced above. 
One observes that the dynamics of the matter variables is given by a time- 
ordered exponential function whose exponent is a bilinear functional of the 
current super-operators J± (x) . Formally we may write Eq. ( p8| ) as 

Pm(t f ) = Tl exp (i£[J+, J_]) p m (ti), (29) 

where we have introduced an influence phase functional 

/■*/ i /■*/ r l s 

i<P[J+,J-}= d i xC c {x) + - d 4 x / dV (30) 

x {-iD F (x - x')„ v J£(x)J$(x') + iD* F (x - x')^Jt{x).r{x') 
+D-(x - x')^J^{x)r{x') + D+(x- 7>)p,J!t(x)JZ(a>)} . 

It should be remarked that the influence phase <P[J+, J-] is both a functional 
of the quantities J± (x) and a super-operator which acts in the space of den- 
sity matrices of the matter degrees of freedom. There are several alternative 
methods which could be used to arrive at an expression of the form (j30j) 
as, for example, path integral techniques [[l] or Schwinger's closed time-path 
method jll). The expression ( |30| ) for the influence phase functional has been 
given in Rcf. [Q without the Coulomb term and for the special case of zero 
temperature. In our derivation we have combined super-operator techniques 
with methods from field theory, which seems to be the most direct way to 
obtain a representation of the reduced density matrix. 

For the study of decoherence phenomena another equivalent formula for 
the influence phase functional will be useful. To this end we define the com- 
mutator function 

D(x-x')^ = i[A^(x),A„{x')} 

= i(D + (i-iV-D_(x-iV) (31) 

and the anti-commutator function 



= D + (x-x')^+D.{x-x')^. (32) 
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in terms of D(x - 




and D\ (x - 










D + {x- 










X )fj,u, 




D_(x- 










X 




iDp{x — 




= \d i{x - 




^sign(< 


-t')D(x- 




—iDp(x — 






)/j.i>' 


-sign(t 


-t')D(x- 





Of course, the previously introduced correlation functions may be expressed 

(33) 
(34) 
(35) 
(36) 

Correspondingly, we define a commutator super-operator J c (x) and an anti- 
commutator super-operator J a {x) by means of 

J£{x) P = \j»(x),p], J»(x)p = {j»(x),p}, (37) 

which are related to the previously introduced super-operators J±(x) by 

J»(x) = J»(x) - Jt(x), J£{x) = Jt(x) + Jl(x). (38) 

In terms of these quantities the influence phase functional may now be written 

as 

i<P[J c ,J a ] = J d^xLdx) (39) 
f d 4 x j dV Ud(x - x')^(x)r a {x') 



- Dl {x-x>)^{x)r c {x>) 



This form of the influence phase functional will be particularly useful later 
on. It represents the influence of the radiation field on the matter dynamics 
in terms of the two fundamental 2-point correlation functions D( x — x') and 
D\(x — a;'). Note that the double space-time integral in Eq. (39) is already 
a time-ordered integral since the integration over t' = x extends over the 
time interval from ti to t = x°. 

For a physical discussion of these results it may be instructive to compare 
Eq. (|2g| ) with the structure of a Markovian quantum master equation in 
Lindblad form ||, 

d Pm -rrr 1 , f A a! 1 A 1 



dt 



i[H m , Pm ] + \Aip m A\ - l -A\A lPm - ~p m A\AiJ , (40) 



where H m generates the coherent evolution and the Ai denote a set of opera- 
tors, the Lindblad operators, labeled by some index i. One observes that the 
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terms of the influence phase functional involving the current super-operators 
in the combinations J + J_ and J_ J + correspond to the gain terms in the 
Lindblad equation having the form Aip m A\. These terms may be interpreted 
as describing the back action on the reduced system of the matter degrees 
of freedom induced by "real" processes in which photons are absorbed or 
emitted. The presence of these terms leads to a transformation of pure states 
into statistical mixtures. Namely, if we disregard the terms containing the 
combinations J+J_ and J_ J + the remaining expression takes the form 



Pm 

(t f ) a U(t f ,ti)pm(U)U*(tf,ti), 



where 

U(t f ,U)=Tiexp 



d 4 xH c (x) 



d A x dVfl f (i-i'W(i)j , '(i') 



(41) 



(42) 



Eq. ( f4l|) shows that the contributions involving the Feynman propagators 
and the combinations J+J+ and J_ J_ of super-operators preserve the purity 
of states jl2| . Recall that all correlations functions have been defined in terms 
of the transversal radiation field. We may turn to the covariant form of the 
correlation functions if we replace at the same time the current density by 
its transversal component j^.. The expression (42) is then seen to contain the 
vacuum-to- vacuum amplitude A[j] of the electromagnetic field in the presence 
of a classical, transversal current density jfc (x) M] , 



A[j] = exp 



-~ J d 4 x J d 4 x'D F (x-x')^j£(x)X T (x') 



(43) 



With the help of the decomposition (|35]) of the Feynman propagator into a 
real and an imaginary part we find 



A[j] = exp [i (S*- 1 



i5« 



(44) 



The vacuum-to- vacuum amplitude is thus represented in terms of a complex 
action functional with the real part 

5 (D = ifd A x f dVsign(i - t')D(x - ^jfiKM, (45) 



and with the imaginary part 

5(2) = \Jd 4 xJ dVDiOr-zV^KOO- 



(46) 



The imaginary part S^ 2 ' yields the probability that no photon is emitted by 
the current j£ t , 



\A[j}\ 2 



exp 



(47) 
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In covariant form we have 

D{x - xV = -7^sign(i - t')S[(x - x') 2 ]g^, (48) 

and, hence, 

S {1) = ~ Jd'xj dVS[(z - x') 2 ]jZ(x)j»(x'). (49) 

This is the classical Feynman- Wheeler action. It describes the classical mo- 
tion of a system of charged particles by means of a non-local action which 
arises after the elimination of the degrees of freedom of the electromagnetic 
radiation field. In the following sections we will demonstrate that it is just 
the imaginary part S^ 2 ' which leads to the destruction of coherence of the 
matter degrees of freedom. 



4 The Interaction of a Single Electron with the 
Radiation Field 

In this section we shall apply the foregoing general theory to the case of a 
single electron interacting with the radiation field where we confine ourselves 
to the non-relativistic approximation. It will be seen that this simple case 
already contains the basic physical mechanism leading to decoherence. 



4.1 Representation of the Electron Density Matrix in the 
Non-Relativistic Approximation 



The starting point will be the representation (29) for the reduced matter 
density with expression (|3^) for the influence phase functional <P. It must be 
remembered that the correlation functions D{x — x')^ l/ and Di{x — x')^ v have 
been defined in terms of the transversal radiation field using Coulomb gauge 
and that they thus involve projections onto the transversal component. In 
fact, we have the replacements, 

D(x - — > D(x - x%j = - ^6 i: j - -^-^ D(x - x') 

for the commutator functions, and 

Di(x - xV — » D x {x - x% = + (dij - Diix - x') 

for the anti-commutator function, where 

/d 3 A; 
2(2^)3^ t cxp (~ ik ( x ~ x ')) ~ ex P ( ik ( x ~ x '))] ' 
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and 

f d 3 k 

Dx{x -x') = J 2(2 ^ [exp {-ik{x - x 1 )) + exp (ik(x - x'))] coth (/3w/2) , 

with the notation — (ui,k) — (|fe|,fe) for the components of the wave 
vector. It should be noted that the commutator function is independent of the 
temperature, while the anti-commutator function does depend on T through 
the factor coth(/3o;/2) = 1 + 2N(u>), where N(u>) is the average number of 
photons in a mode with frequency uj. Hence, invoking the non-relativistic 
(dipole) approximation we may replace 



D(x-x% — > D(t-t% = S i:j D(t-t') = S i:j / duJ(u))smw(t-t'), (50) 

Jo 

and 

D x {x - x')ij — >D 1 (t- t')ij = SijD^t - t') (51) 

/"OO 

— 5ij / dujJ(u) coth {(3lo/2) cos u>(t — <'), 



o 

where we have introduced the spectral density 

J(w) = ^jwe(«-«), (52) 

with some ultraviolet cutoff f2 (see below) . It is important to stress here that 
the spectral density increases with the first power of the frequency u). Had 
we used dipole coupling — ex ■ E of the electron coordinate x to the electric 
field strength E, the corresponding spectral density would be proportional 
to the third power of the frequency. This means that the coupling to the 
radiation field in the dipole approximation may be described as a special 
case of the famous Caldeira-Leggett model g] and that in the language of 
the theory of quantum Brownian motion [Q the radiation field constitutes a 
super-Ohmic environment |L5]Jl(|. Note also that we now include the factor 
e 2 into the definition of the correlation function. Within the non-relativistic 
approximation we may thus replace the current density by 

j(t, x) — > A S ( X - X (t)) + S(x - x(t))^-, (53) 
Ira 2m 

where p(t) and x(t) denote the momentum and position operator of the 
electron in the interaction picture with respect to the Hamiltonian 

Hrn = ^+ V(X) (54) 

2m 



for the electron, V(x) being some external potential. 
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We are thus led to the following non-relativistic approximation of Eq. ( |29| 
Pm(tf) = TV- ( exp 



dt I dt' {-D(t-t') Pc{ t )Pa{t ' ) (55) 
2 mm 



2 mm 



Pm{U). 



This equation represents the density matrix (neglecting the spin degree of 
freedom) for a single electron interacting with the radiation field at tempera- 
ture T. In accordance with the definitions ( |3^ ) and ([38j) p c is a commutator 
super-operator and p a an anti-commutator super-operator. In the theory of 
quantum Brownian motion the function D{t — t') is called the dissipation 
kernel, whereas D\(t — t') is referred to as noise kernel. 

4.2 The Path Integral Representation 



The reduced density matrix given in Eq. (|5^) admits an equivalent path 
integral representation |]l4| which may be written as follows, 

p m (x f ,x'f,tf) = J d 3 Xi J d 3 x' l J(xf 1 x' f ,tf 1 x il x' ll ti)p m (x i ,x' i ,t i ), (56) 

with the propagator function 

J(xf,x' f ,tf,x i ,x' i ,t i ) = J DxDx' exp {i (S m [x] - S m [x'}) + i<P[x,x'}} . 

(57) 

This is a double path integral which is to be extended over all paths xit) and 
x'{t) with the boundary conditions 

x'(t i ) = x' i , x'(tf) = x' f , x(ti) = Xi, x(tf) = x f . (58) 

S m [x] denotes the action functional for the electron, 

S m [x]=f / dt Qmi 2 - V(x)\ , (59) 

while the influence phase functional becomes, 

i#[x, x'] = J ' dt J dt' | iD(t - t') (x(t) - x'(t)) (x(t') + x'{t')) 

—D^t - H) (x(t) x'(t)) (x(t>) - x'(t')) | . (60) 

We define the new variables 

q = x — x', r = —(x + x'), (61) 
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and set, for simplicity, the initial time equal to zero, = 0. We may then 
write Eq. (E§) as 



Pm(rf,qf,tf) = I ' A 3 r, J A 3 q l J(r f ,q f ,t f ;r l ,q. l )p m (r l ,q l ,0). (62) 

The propagator function 

J(rf,qf,tf,ri,qi)= J Br J Dqexp{L4[r, q}} (63) 

is a double path integral over all path r(t), q(t) satisfying the boundary 
conditions, 

r(0)=r h r(t f )=r f , q(0) = q h q(t f ) = q f . (64) 

The weight factor for the paths r(t), q(t) is defined in terms of an effective 
action A functional, 

A[r, q] = j ' At (mrq - V(r + ^q) + V(r - ^q)^j 

+ [ ' At [ 1 At'0{t - t')D{t - t')q{t)r{t') 
Jo Jo 

rt, rt, 



+ - At At , D 1 (t-t')q(t)q(t'). 



(65) 



'o Jo 

The first variation of A is found to be 



5A=- dti Sq(t) 



mr( t ) + -V r (V(r + -q) + V(r - -q)) 



+±J* At'Dlt - t')r{t') + ~f dt>D 1 (t t')q(t>) 



+Sr{t) 



mq{t) + 2V q {V(r + ±q) + V(r - ±q)) 



Aj\t'D(t'-t)q(t')jY 



A_ r*t 

which leads to the classical equations of motion, 



(66) 



mr(t) + \v r {V{r + \q) + V(r - \q)) + A J dt ' D (t - t')r(t') 

= -~j\t l D l {t-t')q{t'), (67) 



and 



iq(t) + 2V q (V(r + ±q) + V(r - ±q)) + ±J * At'D(t' - t)q(f) = 0. (68) 
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4.3 The Abraham-Lorentz Equation 

The real part of the equation of motion ( |67j ) , which is obtained by setting the 
right-hand side equal to zero, yields the famous Abraham-Lorentz equation 
for the electron fll7|. It describes the radiation damping through the damping 
kernel D(t — t') Jl5|. To see this we write the real part of Eq. (p7T) as 



d f* 

mr{t) + — / dt'D(t - t')r(t') = F c 
dt J 



tit), 



(69) 



where F cxt (t) denotes an external force derived from the potential V. The 
damping kernel can be written (see Eqs. (^0|) and (|5^)) 

f n e 2 e 2 d f n 

D(t-t')= dcu — T wsinw((- t') = — ^— dco cos iv(t ~ t') 
Jo ^ 3?H di' Jo 



d amO(t-tf) _ 



d 



3tt 2 dt' t - t' 3tt 2 dt 

where we have introduced the function 

/(*) - 

To be specific the UV-cutoff fi is taken to be 



which implies that 



where 



Tifi = mc , 



ft Ac 



Ac 



ft_ 

mc 



(70) 

(71) 
(72) 
(73) 



is the Compton wavelength. For an electron we have 

A c ~ 3.8 x 10~ 13 m and Q w 0.78 x 10 21 s -1 . (74) 



The term of the equation of motion ( p9[ ) involving the damping kernel 
can be written as follows, 



d e 2 d /"* 

■f / dt'D(t-t')r(t') = ^^- / dt' 
dtj v ' y 1 37r 2 dtJ 



r(f') (75) 



e" a 
3^dt 



dt'f(t-t')f(t') + f(0)r(t)-f(t)r(0). 



For times t such that fit^l, i.e. £ 10 21 s, we may replace 

fit) _^ nS ( t ), 



(76) 
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and approximate f(t) « 0, while Eq. (|7Cj) yields /(0) = fl. Thus we obtain, 

e 2 

I A+'na - t'M-(t') = - 

dt 

which finally leads to the equation of motion, 



j\t'D{t-t')r{t') = ^± [-lr(t) + nr(t)} , (77) 



e ^y m -f- m = F^{tl (78) 

where v = r is the velocity. This is the famous Abraham-Lorentz equation 
| jl7| . The term proportional to the third derivative of r (t) describes the damp- 
ing of the electron motion through the emitted radiation. This term does not 
depend on the cutoff frequency, while the cutoff-dependent term yields a 
renormalization of the electron mass, 

mji = to + Am = m H (79) 

3tt 2 

It is important to note that the electro-magnetic mass Am diverges linearly 
with the cutoff. The equation of motion ( |78| ) can be obtained heuristically by 
means of the Larmor formula for the power radiated by an accelerated charge. 
More rigorously, it has been derived by Abraham and by Lorentz from the 
conservation law for the held momentum, assuming a spherically symmet- 
ric charge distribution and that the momentum is of purely electromagnetic 
origin Jl7[ . 

For the cutoff fl chosen above we get 

tug 2 4 

and, hence, 

Am 4 

= —a w 0.0031. 

m 6tt 

The decomposition ( |79| ) of the mass is, however, unphysical, since the electron 
is never observed without its self-held and the associated held momentum. In 
other words, we have to identify the renormalized mass with the observed 
physical mass which enables us to write Eq. ( |7S| ) as 

m R [«(*) - r v(t)] = F cxt (t). (81) 

Here, the radiation damping term has been written in terms of a characteristic 
radiation time scale To given by 

to = - \r e w 0.6 x 10- 23 s, (82) 

D7T777fl 3 

where r e denotes the classical electron radius, 

a\ c « 2.8 x 10" 15 m. (83) 



47T777 



R 
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It is well-known that Eq. (pl|), being a classical equation of motion for the 
electron, leads to the problem of exponentially increasing runaway solutions. 
Namely, for F cxt = we have 

v - t v = 0. (84) 

In addition to the trivial solution of a constant velocity, v — const, one also 
finds the solution 

v(t) = v(0) exp^/ro), 

describing an exponential growth of the acceleration for i)(Q) ^ 0. In order 
to exclude these solutions one imposes the boundary condition 

v(t) — > for t — > oo, 

if ,F cx t also vanishes in this limit. This boundary condition can be imple- 
mented by rewriting Eq. (|8l]) as an integro-differential equation 

POO 

m R r(t) = / dsexp(-s)F oxt (i + r s). (85) 
Jo 

On differentiating Eq. j85| ) with respect to time, it is easily verified that one 
is led back to Eq. (fn]). However, for -F C xt = it follows immediately from 
Eq. ( |S5"| ) that v = const, such that runaway solutions are excluded. 

On the other hand, Eq. (^H|) shows that the acceleration depends upon 
the future value of the force. Hence, the electron reacts to signals lying a time 
of order tq in the future, which is the phenomenon of pre-acceleration. This 
phenomenon should, however, not be taken too seriously, since the description 
is only classical. The time scale To corresponds to a length scale r e which 
is smaller than the Compton wavelength Ac by a factor of a, such that a 
quantum mechanical treatment of the problem is required. 



4.4 Construction of the Decoherence Function 

In this subsection we derive the explicit form of the propagator function ( p3| ) 
for the reduced electron density matrix in the case of quadratic potentials, 

V(x) = ^m R Lj 2 x 2 . (86) 

Our aim is to introduce and to determine the decoherence function which 
provides a quantitative measure for the degree of decoherence. On using 

V{r + q/2) + V{r - q/2) = m R u; 2 r 2 + m R Lu 2 q 2 /A (87) 
-V(r + q/2) + V(r - q/2) = -m R cu 2 r ■ q, 
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the classical equations of motion take the form 

/>oo 

r(i)+cjo / dsexp(-s)r(t + T s) 
Jo 

/>oo 

q{t) + / ds exp(-s)q(t - t s) 
Jo 



m R 



tor 



= 0. 



— jTdt'ACt-fj^Xss) 

(89) 



Note, that Eq. ( |89| ) is the backward equation of the real part of Eq. (B8|). More 
precisely, if g(£) solves Eq. j8S|), then r(t) = q(tf — t) is a solution of ( p8| ) 
with the right-hand side set equal to zero. 

The above equations of motion lead to the following renormalized action 
functional 



A[r, q] = 



dt m r 



r(t)q(t) - u%q(t) / dsexp(-s)r(i + r s) 



dt / dt'Di(t-t')q{t)q(t'). 



(90) 



In the following we shall use this renormalised action functional instead of 
the action given in Eq. ([35]). By variation with respect to q(t) we immediately 
obtain Eq. (|S8"1), whereas the variation with respect to r(i) yields: 



dt mn 



q(t)Sr(t) + loI 



/•oo 

/ dsexp(— s)q(t)5r(t + tqs) 
■/o 



0, 



which implies 



dtq(t)5r(t) + uo^ 



POO rtf 

I ds I dtcxp(— s)q(t)5r(t + tos) 
Jo Jo 



dtq(t)5r(t) + ul 



tf+T S 



dt exp(—s)q(t — T s)Sr(t) 



= 0. 



(91) 



In the last time integral we may extend the integration over the time interval 
from to tf. This is legitimate since to is the radiation time scale: By setting 
this variation of the action equal to zero we thus neglect times of the order 
of the pre-acceleration time, which directly leads to the equation of motion 

©• 

Since the action functional is quadratic the propagator function can be 
determined exactly by evaluating the action along the classical solution and 
by taking into account Gaussian fluctuations around the classical paths. We 
therefore assume that r(t) and q(t) are solutions of the classical equations 
of motion (^8|) and (|8^) with boundary conditions (Q). The effective action 
along these solutions may be written as 



Alb", 9] = m R [r f q f - 
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dtm R q(t) 



f(t)+uj I dsexp(-s)r(t + tqs) 



4 Jo 



dt / dt'D x {t-t')q(t)q{t'), 



(92) 



or, equivalently, 



Ad[r,q] =m R [r f q f -r-iqi 



dt q(t) 



dt 



dt' 'Dx(t - t')q(t') 



tf 



dt / dt'Di(t-t')q(t)q(t'). 



(93) 



Eq. ( p8[ ) shows that the solution r(t) is, in general, complex due to the cou- 
pling to q(t) via the noise kernel Di(t — t'). Consider the decomposition of 
r(t) into real and imaginary part, 

r(t) =r^(t) +ir^(t), (94) 

where is a solution of the real part of Eq. (^8|), while solves its 
imaginary part, 



m R 



r (2) '(t) + to 2 / dsexp(-s)r (2) (t + t s) 



i d n 

2~dl 



dt'D^t-t^qit'). 

(95) 

We now demonstrate that, in order to determine the action along the 
classical paths, it suffices to find the homogeneous solution and to insert 
it in the action functional 11141 . In other words we have 



where 



A c i[r^,q]=A c i[r,q], 



A d [r^ 1 \q]=m R [rfq i -rf ) q i ] 



(96) 



i 

"4 Jo 



dt 



tf 



To proof this statement we first deduce from Eq 



dt'Dxit-t^qitfqit'). 
b that 



(97) 



i r 1 d r 1 

2 1 dtq{t) dl I dt ' D ^ t - t ')^') 



im R dtq{t) r {2) \t) + u? Q \ ds exp(-s)r (2) (t + t s) 
Jo Jo 



-vm R [ryqf - rf ] qi] 



f l f \ f°° 
-\m R I dt r (2) (t)q{t) + u 2 / ds exp(-s)r (2) (t + T s)q(t) 
Jo I Jo 



.(98) 
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The term within the square brackets is seen to vanish if one employs Eq. ( p9| ) 
and the same arguments that were used to derive the equation of motion 
from the variation of the action functional. Furthermore, we made use 
of r*( 2 )(0) = r( 2 \tf) = which means that the real part r^(t) of the solution 
satisfies the given boundary conditions. Hence we find 



\ £ dt q{t) ^t Si* dt ' Di{t ~ t,)q{tl) = - imR ^? )q f - ^ )q ^ (99) 

from which we finally obtain with the help of j9^), 

Ad[r, q] = m R [r f q f - r.q,} - im R [rf } q f - rf'qi] 

+ - f f dt ( * tib'D^t-t'WW) 
4 Jo Jo 

= Mr {1) ,q]- (100) 

This completes the proof of the above statement. 

Summarizing, the procedure to determine the propagator function for the 
electron can now be given as follows. One first solves the equations of motion 

poo 

r(t)+wg / dsexp(-s)r(t + T s) = 0, (101) 
Jo 

POO 

q{t)+ul \ dsexp(-s)q{t-T s) = 0, (102) 
Jo 

together with the boundary conditions (p4|). With the help of these solutions 
one then evaluates the classical action, 

Ad[r, q] = m R [r f q f - r.q,} + i f ' dt [ ' dt'D x {t - t')q(t)q(t'), (103) 

4 Jo Jo 

which immediately yields the propagator function 

J{rf,qf,tf,ri,qi) = Nexp {iA c i[r, q]} 

= Nexp{im R ,(rfq f - r % q{) + r(q f ,qi,t f )} . (104) 

Here, N is a normalization factor which is determined from the normalization 
condition 

J d 3 r f J{r f ,q f = O.f/jr^gO = 5(«n). (105) 

The function r(q/,qi,tf) introduced in Eq. (104) will be referred to as the 
decoherence function. It is given in terms of the noise kernel D\(t — t') as 

r(q f , q h t f ) = -- T dt T dt'D^t - t')q(t)q(t'). (106) 
4 Jo Jo 
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Explicitly we find with the help of Eq. (|T]) , 



r = -- 



4 Jo 



dt 



dt' 



dujj(cu) coth(/3w/2) cos uj(t - t')q(t)q{t'). (107) 



The double time-integral can be written as 



Re/ dt dt' exp[iu(t - t')]q(t)q(t') = / dt exp(iujt)q(t) 
Jo Jo Jo 

Hence, the decoherence function takes the form 

1 f°° 

r(q f , qi ,t f )=-- J d W JHcoth(/3 W /2)|Q( W )| 2 , 



where we have introduced 

Q(u) 



dt exp(iuit)q(t). 



(108) 



(109) 



(110) 



It can be seen from the above expressions that _T is a non-positive func- 
tion. The decoherence function will be demonstrated below to describe the 
reduction of electron coherence through the influence of the radiation field. 



5 Decoherence Through the Emission of 
Bremsstrahlung 

As an example we shall investigate in this section the most simple case, 
namely that of a free electron coupled to the radiation field. This case is 
of particular interest since it allows an exact analytical determination of 
the decoherence function and already yields a clear physical picture for the 
decoherence mechanism. Having determined the decoherence function, we 
proceed with an investigation of its influence on the propagation of electronic 
wave packets. 



5.1 Determination of the Decoherence Function 



We set ujq = to describe the free electron. The equations of motion (101) 
and (102) with the boundary conditions (|64|) can easily be solved to yield 



r(t) = n + 



Tf-Ti 



q(t) 



If 



tf 



(111) 



Mak ing use of Eq. (104) and determining the normalization factor from 
Eq. ( |105[ ) we thus get the propagator function, 



J(rf,qf,tf,ri,qi) 
( tor \ 3 (.m R 



: (r f - r t ){q f - q % ) + r(q f ,q h t f ) 



(112) 
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As must have been expected J is invariant under space translations since it 
depends only on the difference Vf — rj. Furthermore, one easily recognizes 
that the contribution 

G(rf-ri,q f -qi,tf)= (jj^J exp |i^(r/ - r»)(qr/ - (113) 

is simply the propagator function for the density matrix of a free electron 
with mass ttir for a vanishing coupling to the radiation field. We can thus 
write the electron density matrix as follows, 

Pm{rf,qf,t f ) = J <^ r i J d 3 q l G(r f - n,q f - q i: tj) 

x enp{r(q f ,qi,t f )} p m (ri,qi,0), (114) 

which exhibits that the decoherence function r describes the influence of the 
radiation field on the electron motion. 

We proceed with an explicit calculation of the decoherence function. It 
follows from Eqs. ( |TTo| ) and (JTT^) that 

r>t \ [ tf A+ r j.\Qf-9i exp(iorf/)-l 

Qw) = I dicxp(icjt) — = w, (115) 

Jo tf 1W 

where 

w = —(q f -qi). (116) 
*/ 

Therefore, the decoherence function is found to be 

r = -^£ /"^ i-cos^ coth(/3a)/2)) (U7) 

6tt z J w 

where we have used expression (|5^) for the spectral density J(uS). The deco- 
herence function may be decomposed into a vacuum contribution _T vac and a 
thermal contribution _T t h, 

r = r vac + r th , (us) 

where 

e 2 w 2 f° , 1-coswtf 
r vac = — — =- / dw L (119) 

r th = /" duj 1 —^^- [coth(/3c/2) - 1] . (120) 

6tt z Jq uj 

The frequency integral appearing in the vacuum contribution can be eval- 
uated in the following way. Substituting x — u>tf we get 

/ da; dx =h\ntt + C + 0[-—\, (121) 

Jo UJ Jo x \MtfJ 



and 
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where C ss 0.577 is Euler's constant Jig] , For Qtf 3> 1 we obtain asymptoti- 
cally 



e 2 u> 2 , „. e 2 , „. (q f - qi) 2 



Gtt 



(122) 



To determine the thermal contribution T^h we first write Eq. (12C) as 
follows, 



r, 



e 2 w 2 /■*/ 



6tt 2 



rh 



dui [coth(/3w/2) — 1] sinwi 



2 2 

e ur 



(123) 



Introducing the integration variable x = (3oj we can cast the double integral 
/ into the form 



I 



1 



dt 



HQ 



dx [coth(x/2) - 1] sin (tx/j3) 



Here, we have f3fi = TiD/hsT and, using the cutoff hfi = mc 2 , we get 

™„2 



k B T' 



For temperatures T obeying 



ksT <C mc 



(124) 



the upper limit of the x-integral may be shifted from (3D to oo. Condition 
( |124|) states that 

h 2 h 2 

> 



mksT m 2 c 2 ' 



which means that the thermal wavelength A t h = h/\/2mkBT is much larger 
than the Compton wavelength, 



Ath ^> Ac- 



(125) 



Thermal and Compton wave leng th arc of equal size at a temperature of 
about 10 9 Kelvin. Condition (|l24| ) therefore means that T <C 10 9 K. Under 
this condition we now obtain 



1 

PJo 
1 r*t 



dt 



dt 



PJo 

fsuMntf/m 



da; [coth(x/2) - 1] sin (te//3) 

7TCOth 



(126) 
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where we have employed the formula 



f°° 1 

/ da; fcoth(x/2) — 1] sin tx = n coth(7rr) . (127) 

Jo t 



The quantity 



7T 



2.4- 1CT 12 s/T[K] 



(128) 



irksT 

represents the correlation time of the thermal radiation field. Putting these 
results together we get the following expression for the thermal contribution 
to the decoherence function, 



th 



6tt 2 



■In 



smh(t f /T B )\ (qf-qi) 

tf/TB 



t 1 



(129) 



Adding this expression to the vacuum contribution (122) and introducing 
a = e 2 /Airhc and further factors of c, we can finally write the expression for 
the decoherence function as 



2a 

r(qf,qi,t f ) fa - — 



In Qtf + In 
Alternatively, we may write 



r(qf,qi,t f ) = 
where the quantity L(tf) defined by 



3?r 



4a 



In Qtf + hi 



sinh(f//r 5 ) 

tf/TB 



(g/ - q*) 2 

2L(t f f ' 



smh(tf/T B ) 

tf/TB 



(Of - 9») 2 
(ctf) 2 



{ctff 



(130) 



(131) 



(132) 



may be interpreted as a time-dependent coherence length. 

The vacuum contribution _T vac to the decoherence function (130) appar- 



ently diverges with the logarithm of the cutoff Q. This is, however, an artificial 
divergence which can be seen as follows. The decoherence fu nctio n is d efine d 
in terms of the Fourier transform Q(u>) of q(t), see Eqs. (109) and (110). 



Evaluating Q(u>) as in Eq. (115) we assume that the velocity is zero prior to 
the initial time t = 0, that it suddenly jumps to the value given by Eq. (116), 
and that it again jumps to zero at time tf. This implies a force having the 
shape of two S- function pulses around t — and t = tf. Such a force acts 
over two infinitely small time intervals and leads to sharp edges in the clas- 
sical path. More realistically one has to consider a finite time scale t v for the 
action of the force which must be still large compared to the radiation time 
scale To- We may interpret the time scale t p as a preparation time since it 
represents the time required to prepare the initial state of a moving electron. 
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A natural, physical cutoff frequency of the order fl ~ l/r p is thus introduced 
by the preparation time scale t p and we may set 

Qt f = ^ (133) 

Tp 

in the following. It should be noted that the weak logarithmic dependence 
on Q shows that the precise value of the preparation time scale t p is rather 
irrelevant. The important point is that the preparation time introduces a 
new time scale which removes the dependence on the cutoff. The vacuum 
decoherence function can thus be written, 

showing that it vanishes for large times essentially as tj 2 . 

The thermal contribution r t h is determined by the thermal correlation 
time tb- For T — ► we have tb — ► oo, and this contribution vanishes. For 
large times tf 3> tb the thermal decoherence function may be approximated 

by 

A h «-^%^, (135) 
Sir tb (ctf) 

which shows that i~th vanishes as tj . Thus, for short times the vacuum 
contribution dominates, whereas the thermal contribution is dominant for 
large times. Both contributions F vac and Jth are plotted separately in Fig. 1 
which clearly shows the crossover between the two regions of time. 



Eq. (132) implies that the vacuum coherence length is roughly of the order 

L(t f ) vac ~c-tf. (136) 

To see this let us assume a typical preparation time scale of the order t p ~ 
10 -21 s. If we take tf to be of the order of Is we find that ln(t//r p ) ~ 48. 
In the rather extreme case tf ~ 10 17 , which is of the order of the age of the 



universe, we get ln(t//r p ) ~ 87. On using 37r/4a w 322 and Eq. (132) for 



T = one is led to the estimate (136) 



5.2 Wave Packet Propagation 

Having obtained an expression for the decoherence function r we now proceed 
with a detailed discussion of its physical significance. For this purpose it will 
be helpful to investigate first how r affects the time-evolution of an electronic 
wave packet. We consider the initial wave function at time t = 0, 



( 1 



3/4 



exp 



(sc - af 



iko(x — a) 



(137) 



26 



Heinz-Peter Breuer and Francesco Petruccione 



x 10 
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100 



Fig. 1. The vacuum contrib utio n _T vac and the thermal contribution i~th of the de- 
coherence function r (Eq. (130)). For a fixed value \qj — qi\ = 0.1 ■ ctb, the two 



contributions are plotted against the time t / which is measured in units of the ther- 
mal correlation time tb- The temperature was chosen to be T = IK. One observes 
the decrease of both contributions for increasing time, demonstrating the vanishing 
of decoherence effects for long times. The thermal contribution Zth vanishes as tj 1 , 
while the vacuum contribution _T vac decays essentially as tj 2 , leading to a crossover 
between two regimes dominated by the vacuum and by the thermal contribution, 
respectively. 



describing a Gaussian wave packet centered at x = a with width cto- With 
the help of Eqs. flll3| ), (114) and (131) we get the position space probability 
density at the final time t / , 



p m (r f ,t f ) ee p m {r f ,q f = 0,t f ) 



{r f - n)qi 



2L{t f f 



The Gaussian integrals may easily be evaluated with the result, 



(138) 



Pm(r f ,t f ) 



3/2 

, , exp 



fry - bf 

2a(t f y 



(139) 
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where 

b^a-^l (140) 
m R 

and 

^/) 2 ^o + i J^ + 7^%- (141) 

This shows that the wave packet propagates very much like that of a free 
Schrodinger particle with physical mass m R . The centre b of the proba- 
bility density moves with velocity —ko/mn, while its spreading, given by 



Eq. (141), is similar to the spreading cr(tf) 2 ICC which is obtained from the 



free Schrodinger equation, 

t 2 

*(*/)Le = 4 + T-hl- ( 142 ) 



Am 2 R a 2 ' 



If we write 



t) f Aa 2 
Am 2 R a 2 [ + L(t f y 



we observe that the decoherence function affects the probability density only 
though the widt h a( tf) and leads to an increase of t he spreading. In view 



of the estimate ( 136 ) the correction term in Eq. (143) is, however, small for 
times satisfying 

L(t f ) ~c-t/ >cr . (144) 

This means that the influence of the radiation field can safely be neglected 
for times which are large compared to the time it takes a light signal to travel 
the width of the wave packet. 



2a 



Fig. 2. Sketch of the interference experiment used to determine the decoherence 
factor. Two Gaussian wave packets with initial separation 2a approach each other 
with opposite velocities of equal magnitude v — ko/mn. 



Let us now study the evolution of a superposition of two Gaussian wave 
packets separated by a distance 2a. This case has been studied already by 
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Barone and Caldeira Jig ] who find, however, a different result. We assume 
that the packets have equal widths oo and that they are centered initially at 
x = ±a = ±(a, 0, 0). The packets are supposed to approach each other with 
the speed v — ko/rriR > (see Fig. 2). For simplicity the motion is assumed 
to occur along the ir-axis. Thus we have the initial state 



ip (x) = A x 



Ad 





3/4 


j exp 


(hi 




3/4 


j exp 





(x — a) 
(x + a) 



lk^ix — a) 



\ko(x 



(145) 



where fco = (fco, 0,0) and A±, A2 are complex amplitudes. Our aim is to 
determine the interference pattern that arises in the moment of collision of 



the two packets at x = 0. Using again Eqs. (113), (114) and (131) and doing 
the Gaussian integrals we find 



Pm(r f ,t f ) 



1 



3/2 

, , exp 



2*v(t f )*J ~ e [ 2a(tf) 2 _ 
{lAJ 2 + \A 2 \ 2 + 2ReA 1 A*ex V [<p{r f )]} . (146) 



We recognize a Gaussian envelope centered at rj — with width cr(tf), an 
incoherent sum \Ai\ 2 + 1 ^4_2 1 2 , and an interference term proportional to AiA^. 
The interference term involves a complex phase given by 



<f( r f) = -2iko»v(l - e) 



2a 2 



r(l-e)- 



(147) 



The term — 2ifco r / describes the usual interference pattern as it occurs for a 
free Schrodinger particle, while the contribution 2ifco r /£ leads to a modifica- 
tion of the period of the pattern. The final time tf = amn/ko is the collision 
time and 

a fco 
tf m R 

is the speed of the wave packets. The factor e is given by 



(148) 



m%L{t f Ya{hf 



4a 2 



f2 



(149) 



Obviously w e alw ays have < s < 1. Furthermore, for the situation consid- 
ered in Eq. (144) we have e<l. Thus, we get 



2a 2 



L{tff 



(150) 
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The last expression clearly reveals that the real part of the phase <p(rf) 
describes decoherence, namely a reduction of the interference contrast de- 
scribed by the factor 





f (2a) 2 " 




exp 




= exp 



distance 



2 (coherence length) 2 



(151) 



which multiplies the interference term. As was to be expected from the general 
formula for r, the decoherence factor D is determined by the ratio of the 
distance of the two wave packets to the coherence length. 

Alternatively, we can write the decoherence factor in terms of the velocity 



(148) of the wave packets. In the vacuum case we then get 

J7T V 7V, / V cJ 



D vac = exp 



(152) 



This clearly demonstrates that it is the motion of the wave packets which 
is responsible for the reduction of of the interference contrast: If one sets 
into relative motion the two components of the superposition in order to 
check locally their capability to interfere, a decoherence effect is caused by 
the creation of a radiation field. As can be seen from Eq. (117) the spectrum 
of the radiation field emitted through the moving charge is proportional to 
1/lo which is a typical signature for the emission of bremsstrahlung. Thus 
we observe that the physical origin for the loss of coherence described by the 
decoherence function is the creation of bremsstrahlung. 



It is important to recognize that the frequency integral of Eq. (117) con- 
verges for uj — ► 0, see Eq. (121). The decoherence function r is thus infrared 
convergent which is obviously due to the fact that we consider here a pro- 
cess on a finite time scale i/. This means that we have a natural infrared 
cutoff of the order of J? m i n ~ V*/' m addition to the natural ultraviolet 
cutoff fi ~ 1/tj, introduced earlier. The important conclusion is that the 
decoherence function is therefore infrared as well as ultraviolet convergent. 

It might be instructive, finally, to compare our results with the corre- 
sponding expressions which are derived from the famous Caldeira-Leggett 
master equation in the high-temperature limit (see, e.g. ||). From the latter 
one finds the following expression for the coherence length 



L{t f ) 



2 _ A th 
CL 



27*/' 



(153) 



where 7 is the relaxation rate. This is to be compared with the expressions 
(132) for the coherence length. For large temperatures we have the following 
dominant time and temperature dependence, 



L{tff ~ g 



T 



and 



(154) 
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Hence, while both expressions for the coherence length are proportional to 
the inverse temperature, the time dependence is completely different. Namely, 
for tf — > oo we have 

L (t f ) 2 — ► oo and L(t f ) 2 ch — > 0, (155) 

and, therefore, complete coherence in the case of bremsstrahlung and total 
destruction of coherence in the Caldeira-Leggett case. 



6 The Harmonically Bound Electron in the Radiation 
Field 



As a further illustration let us investigate briefly the case of an electron in the 
radiation field moving in a harmonic external potential. Another approach 
to this problem may be found in fll9| , where the authors arrive, however, at 
the conclusion that there is no decoherence effect in the vacuum case. 



We take ujo > and solve the equation of motion (101) with the help of 
the ansatz 

r(t) = r cxp(zt), (156) 
where, for simplicity, we consider the motion to be one-dimensional. Substi- 



tuting this ansatz into (101) one is led to a cubic equation for z, 

z 2 -t q z 3 + w 2 = 0. (157) 

For vanishing coupling to the radiation field (tq = 0) the solutions are lo- 
cated at z± = ±iwo, describing the free motion of a harmonic oscillator with 
frequency uq. 

For To > the cubic equation has three roots, one is real and the other 
two are complex conjugated to each other. The real root corresponds to the 
runaway solution and must be discarded. Let us assume that the period of 
the oscillator is large compared to the radiation time, 

r « — . (158) 

U} 

Because of t ~ 10 _24 s this assumption is well satisfied even in the regime 
of optical frequencies. We may thus determine the complex roots to lowest 
order in ujqTo, 

z± = ±iw - -t w 2 . (159) 

The purely imaginary roots ±iwo of the undisturbed harmonic oscillator are 
thus shifted into the negative half plane under the influence of the radiation 
field. The negative real part describes the radiative damping. In fact, we see 
that r(t) decays as exp(— "ft/2), where 

7 = T c 2 = \aM_ (160) 
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is the damping constant for radiation damping J17[ . In the following we con- 
sider times tf of the order of magnitude of one period to^tf ~ 1. Because of 
jtf = (ujoTo)(u)Qtf) we then have jtf ~ tqloq "C 1. In this case the damping 
can be neglected and we may use the free solution in order to determine the 
decoherence function. 

Let us consider again the case of a superposition of two Gaussian wave 
packets in the harmonic potential. The packets are initially separated by a 
distance 2a and approach each other with opposite velocities of equal mag- 
nitude such that they collide after a quarter of a period, tf = tt/2loq. The 
corresponding free solution q(t) is therefore given by 



q(t) = qi cos wot + qf smujot. 



(161) 



To describe the situation we have in mind we take qi — 2a (initial separation 
of the wave packets) and g/ = (to get the probability density). Hence, we 
have 

q(t) = — 2aujQ sinujQt, (162) 
and we evaluate the Fourier transform, 



Q(uj) = / dt exp(iwt)q(t) 
Jo 

exp(i[w + uJo]tf) — 1 exp(i[w — oj ]tf) — 1 



UJ + U) 

This yields the decoherence function 



LU — UJq 



r = r( qf = o,q h t f ) 



(163) 



1 — cos(w 



too) 



uJo)tf 1 — cos(o; 



UJ )tf 



{lo - w ) 2 



coth | — 
2 



We discuss the case of zero temperature. The frequency integral in Eq. 



(163) then approaches asymptotically the value 21nJ7i/ which leads to the 
following expression for the decoherence factor, 



exp r v 



exp 



3tt 



In 



(164) 



The interesting point to note here is that this equation is the same as Eq. (152) 
for the free electron, with the only difference that the square (v/c) 2 of the 
velocity, which was constant in the previous case, must now be replaced with 
its time averaged value ((v/c) 2 ). 



7 Destruction of Coherence of Many-Particle States 



For a single electron the vacuum decoherence factor ( |152[ ) turns out to be very 
close to 1, as can be illustrated by means of the following numerical example. 
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We take r p to be of the order of 10 _21 s and tt of the order of Is. Using a 
velocity v which is already as large as 1/10 of the speed of light, one finds 
that -T V ac ~ 10~ 2 , corresponding to a reduction of the interference contrast of 
about 1%. This demonstrates that the electromagnetic field vacuum is quite 
ineffective in destroying the coherence of single electrons. 

For a superposition of many-particle states the above picture can lead, 
however, to a dramatic increase of the decoherence effect. Consider the su- 
perposition 

\1>) = IVi) + IV*) (165) 

of two well-localized, spatially separated TV-particle states and \4>2)- We 
have seen that decoherence results from the imaginary part of the influence 
phase functional ^[J c , J a ], that is from the last term on the right-hand side 
of Eq. ( |39| ) involving the anti-commutator function D\{x — x')^ iy of the elec- 
tromagnetic field. Thus, it is the functional 

r[J c ] =-J| ! 'd 4 x|''d 4 x'D 1 (i- I ') F JfWJ;(x'), (166) 

which is responsible for decoherence. This shows that the decoherence func- 
tion for A^-electron states scales with the square N 2 of the particle number. 



Thus we conclude that for the case of the superposition (165) the decoherence 
function must be multiplied by a factor of ./V 2 , that is the decoherence factor 
for iV-particle states takes the form, 



~ ex P 



8a, ftf \ /v\ 2 t9 



07T V 7V, / V C / 



(167) 



This scaling with the particle number obviously leads to a dramatic increase 
of decoherence for the superposition of TV-particle states. To give an example 
we take N = 6 TO 23 , corresponding to 1 mol, and ask for the maximal velocity 



v leading to a 1% suppression of interference. With the help of (167) we find 
that v r~j 10~ 16 m/s. This means that, in order to perform an interference 
experiment with 1 mol electrons with only 1% decoherence, a velocity of at 
most 10 _16 m/s may be used. For a distance of lm this implies, for example, 
that the experiment would take 3 x 10 8 years! 



8 Conclusions 

In this paper the equations governing a basic decoherence mechanism oc- 
curring in QED have been developed, namely the suppression of coherence 
through the emission of bremsstrahlung. The latter is created whenever two 
spatially separated wave packets of a coherent superposition are moved to one 
place, which is indispensable if one intends to check locally their capability to 
interfere. We have seen that the decoherence effect through the electromag- 
netic radiation field is extremely small for single, non-relativistic electrons. 
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The decoherence mechanism is thus very ineffective on the Compton length 
scale. An important conclusion is that decoherence does not lead to a local- 
ization of the particle on arbitrarily small length scales and that no problems 
with associated UV-divergences arise here. 

The decoherence mechanism through bremsstrahlung exhibits a highly 
non-Markovian character. As a result the usual picture of decoherence as a 
decay of the off-diagonals in the reduced density matrix does not apply. In 
fact, consider a superposition of two wave packets with zero velocity. The 



expression (131) for the decoherence function together with the estimate 
L(tf) vac ~ c- tf for the vacuum coherence length L(i/) va c show that de- 
coherence effects are negligible for times t / which are large in comparison to 
the time it takes light to travel the distance between the wave packets. The 
off-diagonal terms of the reduced density matrix for the electron do therefore 
not decay at all, which shows the profound difference between the decoher- 
ence mechanism through bremsstrahlung and other decoherence mechanisms 
(see, e. g. @). 

A result of particular interest from a fundamental point of view is that 
coherence can already be destroyed by the presence of the electromagnetic 
field vacuum if superpositions of many-particle states are considered. An im- 
portant conclusion which can be drawn from this picture of decoherence in 
QED refers to various alternative approaches to decoherence and the closely 
related measurement problem of quantum mechanics: In recent years sev- 
eral attempts have been made to modify the Schrodinger equation by the 
addition of stochastic terms with the aim to explain the non-existence of 
macroscopic superpositions through some kind of macrorealism. Namely, the 
random terms in the Schrodinger equation lead to a spontaneous destruc- 
tion of superpositions in such a way that macroscopic objects are practically 
always in definite localized states. Such approaches obviously require the in- 
troduction of previously unknown physical constants. In the stochastic theory 
of Ghirardi, Pearle and Rimini ^T|, for example, a single particle microscopic 
jump rate of about 10 _16 s _1 has to be introduced such that decoherence is 
extremely weak for single particles but acts sufficiently strong for many par- 
ticle assemblies. It is interesting to observe that the decoherence effect caused 
by the presence of the quantum field vacuum yields a similar time scale in a 
completely natural way without the introduction of new physical parameters. 
Thus, QED indeed provides a consistent picture of decoherence and it seems 
unnecessary to propose new ad hoc theories for this purpose. 

It must be emphasized that the above picture of decoherence in QED 
has been derived from the well-established basic postulates of quantum me- 
chanics and quantum field theory. It therefore does not, of course, constitute 
a logical disprove of alternative approaches. However, it does represent an 
example for a basic decoherence mechanism in a microscopic quantum field 
theory. In particular, it provides a unified explanation of decoherence which 
does not suffer from problems with renormalization (as they occur, e.g. in 
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alternative theories [£2| ) and which does not exclude a priori the existence of 
macroscopic quantum coherence. Only under certain well-defined conditions 
regarding time scales, relative velocities and the structure of the state vector, 
it is true that decoherence becomes important. Thus, decoherence is traced 
back to a dynamical effect and not to a modification of the basic principles 
of quantum mechanics. 

In this paper we have discussed in detail only the non-relativistic approxi- 
mation of the reduced electron dynamics. For a treatment of the full relativis- 
tic theory, including a Lorentz invariant characterization of the decoherence 
induced by the vacuum field, one can start from the formal development given 
in section 2. An investigation along these lines could also be of great interest 
for the study of measurement processes in the relativistic domain |p3l . 
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